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The effect of impurities on Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states in a layered super- 
conductor with d-wave pairing symmetry is investigated using the tight-binding model and the 
Bogoliubov-de Gennes equations. At low temperature and a strong exchange or Zeeman field ap- 
plied parallel to its conducting plane, a two-dimensional (2D) square lattice-like Larkin-Ovchinnikov 
state is more energetically favorable in a clean system. In the presence of impurities, the spatial 
profile of the order parameter remains as a 2D square lattice, and it is distorted only near the im- 
purities when the impurity concentration is low. As impurity concentration is increased to a certain 
level, quasi-one-dimensional (ID) like FFLO state becomes more energetically favorable. Increasing 
temperature with fixed impurity concentration can also induce a 2D to ID FFLO state transition. 
Within the present finite size calculation, we did not find the existence of the Fulde and Ferrell state 
before the system becoming normal as the impurity concentration is raised. 



PACS numbers: 74.81.-g, 74.25. Ha, 74.50.+r 

Recently, the possible existence of the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state in heavy fermion com- 
pound CeCoIng^ has revived the interest in this inhomo- 
geneous superconducting state. The FFLO state was pre- 
dicted forty years ago by Fulde and Ferrell 2 , and Larkin 
and Ovchinnikov 3 . for a superconductor subject to a 
strong exchange field or a strong magnetic (Zeeman) field 
applied in the planes of a quasi-two-dimensional super- 
conductor. Soon after the prediction, it was shown that 
the existence of the FFLO state depends on the magni- 
tude of the exchange or Zeeman field and is very sensitive 
to the presence of impurities^. This could be one of the 
reasons that the FFLO state has never been observed in 
conventional low-temperature superconductors. 

CeCoIns is likely to be a quasi-two-dimensional d-wave 
superconductor—. It has been predicted that at low tem- 
perature in a clean two dimensional d-wave supercon- 
ductor, the order parameter of the energetically favored 
Larkin-Ovchinnikov (LO) state has 2D square-lattice-like 
structure&i^, while the ID stripe-like LO state is more 
favored in s-wave superconductors. In the LO state, the 
order parameter is real and spatially inhomogeneous with 
positive and negative signs spreading over the configura- 
tion space. Whereas in the FF state, the order parame- 
ter gets a spatially inhomogeneous phase similar to that 
of a plane-wave, but its magnitude remains to be spa- 
tially uniform. In such a state, a spin current has been 
shown to exist. In the clean limit, CeCoIns is a material 
with very long mean free path. However, it is possible 
to introduce impurities in this material by substituting 
some Ce atoms by La atoms^. The scattering from im- 
purities may break the symmetry of the order parameter 
in the FFLO state. Using Born approximation to treat 
impurity scattering, Agterberg and Yang have derived a 
Ginzburg-Landau (GL) free energy functional for both 
s-wave and d-wave superconductors with non-magnetic 
impurities and Zeeman fields^. Using this free energy, 



they have shown that impurities could induce a change 
in the structure of the FFLO state in a d-wave super- 
conductor, namely, a phase transition from an LO state 
in the clean limit to a FF state with relatively high con- 
centration of impurities. This GL free energy is valid 
near normal to superconductor phase transition line in 
the (T, H) plane. The impurity induced FF phase was 
found to be between the LO phase and the uniform su- 
perconducting phase and thus is away from the normal to 
superconductor phase transition line. Therefore, it is of 
interest to study such a problem with an approach which 
goes beyond the GL theory, and to take the effect due to 
multiple scatterings from impurities into account when 
the impurity concentration is relatively high . 

In this paper, based upon a tight-binding model for a 
d-wave superconductor with FFLO state, we show that 
when impurity concentration is increased to a certain 
level, the spatial profile of the order parameter changes 
from 2D like to ID like. This ID LO phase lies between 
2D FFLO (or LO) phase and the normal phase. No FF 
state before the system becoming the normal state is 
found. In the ID LO phase, the impurities would like to 
reside on or very close to the nodal lines in real space. As 
a result, the nodal lines would not be as straight as those 
in the clean limit. When temperature is raised and with 
fixed impurity concentration, the mean-free path should 
decrease and a 2D to ID LO state phase transition is 
shown to occur as well. 

To investigate the effect of impurities on the FFLO 
state in superconductors, we study a phenomenological 
model on a square lattice: 
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where iy is the hopping integral and /x is the chemi- 
cal potential. Throughout our study, we take to be 
unity for nearest neighbors and zero otherwise. The Zee- 
man energy term ah, with a = ±1 for spin-up and - 
down electrons, arises from the interaction between the 
magnetic field and the spins of electrons. U$ denotes 
the impurity potential. j m , with m = 1 to n, where 
n is the number of impurities. The singlet supercon- 
ducting order parameter has the following definition: 
A i3 =V< c^Cji - c^Cjf > /2. 

This Hamiltonian can be diagonalized by solving the 
discrete Bogoliubov-de Gennes equations: 

where = —Uj - (m + <rh)6ij + U 5i t j m . u™ CT and 

are the Bogoliubov quasiparticle amplitudes on the j-th 
site. The self-consistency condition for the OP: 
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is solved by iteration. Here 7 = (±1, 0) and (0, ±1), and 
At = (A i+£ + Ai_ 2 - A i+ y - A,_y)/4 is the d-SC OP at 
site i. 

We choose Uq = 5.0 for a relatively strong impurity po- 
tential. If the impurity potential is taken to be too high, 
the LO state may be destroyed before we can see a 2D to 
ID transition. The size of the lattice for performing the 
numerical computation is chosen to be 32 x 32. There 
are totally N = 1024 sites. The impurities are randomly 
distributed. For each impurity concentration, we calcu- 
late the order parameter for several independent set of 
impurity configurations. The LO phase transition point 
for a fixed impurity concentration should be obtained by 
taking the average of all trial transition points in differ- 
ent impurity configurations. In order to check whether 
the FF state becomes the ground state when impurity 
concentration is increased, we assign a random phase to 
the order parameter at each site as an initial condition. 

Impurities suppress the superconductivity in a d-wave 
superconductor and they may also change the periodic- 
ity of the order parameter in the LO state. We take 
V = 2.0, /i = —1.0 and h = 0.40 in all the calculations. 
With this set of parameter, the order parameter shows 
2D square-lattice variation for a clean d-wave supercon- 
ductor at low temperature. In our numerical calculation 
with impurities, we assume that the periodicity of the 
LO state along the x direction could be different from 
that of the y direction. 

Configuration space phase transition induced by impu- 
rities in the low temperature limit: Fig. [1] shows the cal- 
culated order parameter for various impurity concentra- 
tions. When the concentration of impurities is very low 
(below n/N w 3.0%), the order parameter still iterates 
practically to that of a 2D lattice LO state as that in a 
clean limit with nodal lines along (110) directions, except 



the distortion occurs near the impurity site, as shown in 
Fig. [I] (a) and (b) for n/N w 1.0% and 2.0%. Impurities 
are found likely to stay on or near nodal lines of the order 
parameter in the LO state. Looking at the spatial profiles 
of the order parameter, one can hardly find the impurity 
sites in these graphs because they stay at places where 
the order parameter is very small or zero. As the impu- 
rity concentration increases, the magnitude of the order 
parameter decreases and at the same time, the distortions 
become stronger and stronger. When the concentration 
reaches a certain level (above 3.5% in our case) as shown 
in Fig. [T] (d) and (e) for n/N w 4.0%, the order pa- 
rameter becomes quasi-one-dimensional. The nodal lines 
no longer cross each other and are approximately along 
(100) direction. Further increasing of impurity concen- 
tration leads to the reduction of the superconducting or- 
der parameter until it vanishes. In the above graphs, the 
impurity site is marked by a "+" in the corresponding 
contour plots. The configurations for the positions of im- 
purities are generated randomly by computers. We also 
have performed the same calculation for 5 other configu- 
rations with the same number of impurities. The essen- 
tial feature of the above conclusion remains practically 
unchanged. For the purpose of illustration, we present 
the results for two different impurity configurations with 
40 impurities in our lattice as shown in Fig. [1] (d) and 
1(e). However, unlike the results in Ref. 10, we did not 
obtain an FF phase in our calculation. The phases of the 
order parameter at different sites are always converging 
to the same value in the self-consistent calculation. We 
have also investigated this effect for an s-wave supercon- 
ductor. In a clean system, the energetically favored state 
is ID stripe-like LO state. Upon increasing impurity con- 
centration, the order parameter remains ID stripe-like 
until it vanishes and becomes normal state. An FF phase 
has never shown up in the present approach. 

Configuration space phase transition in the presence of 
impurities due to the raise of temperatures: To study this 
effect, we fix the number of impurities in the lattice to 
30. The results are shown in Fig. [5] At low temperature 
(T=0.001), the order parameter is clearly having a 2D 
structure even though it is strongly distorted at or near 
the sites of impurities. As the temperature is raised, 
we see the variation becomes more and more ID like. 
At T=0.021, the LO state is destroyed, and the system 
becomes normal. For lower impurity concentration, say 
2.0%, we found no 2D to ID transition. The 2D LO 
state persists up close to the superconducting transition 
temperature where the order parameter vanishes. 

The phase change caused by raising the temperature 
indicates that in a T-H phase diagram, the ID LO phase 
locates between the 2D LO phase and the normal state, 
while the FF state predicted in Ref. 10 lies between the 
LO state and the uniform superconducting state. The 
two studies investigate two different regions in the phase 
diagram. Our result is in the high field region result- 
ing from the fact that our method is more suitable for 
small period FFLO states in strong exchange or Zeeman 
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FIG. 1: Spatial profiles of the order parameter for various 
impurity concentrations. The size of the lattice is 32x32. 
(a) 10 impurities, (b) 20 impurities, (c) 30 impurities, (d) 
40 impurities, and (e) 40 impurities. Here (d) and (e) have 
different impurity configurations. When there are more than 
35 impurities, the nodal lines no longer cross each other. The 
sites of impurities are marked by "+" on the contour plots. 



field, while the study in Ref. 10 might be in relatively low 
field. For low field, the periodicity of the LO state will 
be large, the required lattice size is beyond our numerical 
capability. 

The Local density of states for spin up quasiparticles 
(LDOS) at a site next to an impurity when is impurity 
concentration low: With impurity concentration is in- 




FIG. 2: Spatial variations of Aj at various temperature with 
30 impurities on a lattice of size 32 x 32:(a) T=0.001; (b) 
T=0.005; (c) T=0.01; (d) T=0.02. The impurities are marked 
by "+" on the contour plots. 



creased to 3.5%, our result shows a 2D LO state to ID 
LO state transition. We would like to check whether 
this transition could be observed by the scanning tunnel- 
ing microscopy experiments, which directly measure the 
LDOS of the system and its spatial maps. In Ref. 8, the 
present authors have calculated the LDOS at a site near- 
est neighboring to a single impurity when the impurity is 
located at an saddle point or extreme of the order param- 
eter. The LDOS spectrum at site nearest neighboring to 
the impurity show quite different behavior from a pure 
d-wave superconductor. If the impurity is at an order pa- 
rameter saddle point, the LDOS is practically unaffected 
by the impurity, with no impurity-induced resonant peak 
or peaks appearing. This is because the order parameter 
seen by a quasiparticle before and after scattering by an 
impurity does not change sign, thus there are no zero en- 
ergy Andreev's bound states. If the impurity is at an or- 
der parameter extremum, a ±eo+ (Zeeman energy) pair 
of finite-energy resonance peaks are induced in the LDOS 
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FIG. 3: the LDOS at sites next to three impurity sites when 
there are 10 impurities in a lattice ol size 32x32. V = 2.0, 
h = 0.4. The outside two peaks for each curve, marked by 
red arrows, are coherence peaks. The subgap LDOS are com- 
plicated due to the interferences among the resonance states 
formed in the vicinity of different impurity sites. 

by the impurity, instead of one at near zero energy, as in 
a uniform d-wave superconductor. When multiple im- 
purities are considered, the impurities can not all locate 
at a saddle point or a order parameter extremum. Un- 
less a impurity is located at one of the saddle points, we 
would expect that the LDOS spectrum would be similar 
to the case when the impurity is at an order parameter 
extreme because one can always find a quasiclassical or- 
bit of incident and scattered quasiparticles experiencing 
sign change of the order parameter. 

In fig. [3l we show the LDOS at sites nearest neighbor- 
ing to some impurity sites when impurity concentration 
is low (10 impurities in a lattice of the size 32x32 or the 
impurity concentration is approximately 1%). The sub- 
gap LDOS at these sites are very complicated comparing 
to the case of single impurity. This is due to the interfer- 
ences among the resonance states formed in the vicinity 
of different impurities sites. Again the stable configura- 
tion here for the order parameter is still two-dimensional, 



and the image maps of the LDOS at certain fixed energies 
still exhibit the same checkerboard patterns^. However, 
near the 2D to ID transitional region, the impurity con- 
centration becomes larger and it is between 3 to 4%. Un- 
der this condition, there are strong interferences among 
the states from different impurities and nodal lines. As 
a result, no clear image maps could be obtained to show 
the transition. This conclusion also hold true if we fix 
the impurity concentration at 3%, then raise the tem- 
perature to detect the 2D to ID transition. Since the 
magnetization along the nodal lines should be larger—, 
the 2D to ID transition could be observed by neutron 
scattering experiments. 

In summary, the effect of non-magnetic impurities on 
the FFLO state in a d-wave superconductor is studied 
by using the tight-binding model and the Bogoliubov-de 
Gennes equations. We have shown that the impurities 
can induce a 2D to ID configuration space phase transi- 
tion in such a system. At low temperature, with a strong 
magnetic field applied to its conducting plane, a 2D lat- 
tice like Larkin-Ovchinnikov state is more energetically 
favorable in a clean system. When impurity concentra- 
tion is very low, the spatial profile of the order parameter 
remains practically unchanged, except it is distorted near 
the impurity sites. As impurity concentration is raised to 
a certain level, the scatterings from impurities break the 
symmetry of the system and ID-stripe like LO state be- 
comes more energetically favorable. The ID-stripes are 
not straight because they have to accommodate the im- 
purities to reside on or near the nodal lines. The 2D to 
ID LO transition is also obtained if the impurity con- 
centration is fixed at, say 3.0% and the temperature 
is raised. In the present self-consistent calculations on 
square lattice with size 32x32, we are not able to find the 
FF state as predicted in Ref. 10. 
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